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We consider a dynamical system which has a stable attractor and which is perturbed by an additive 
noise. Under some quite typical conditions, the fluctuations from the attractor are intermittent and 
have a probability distribution with power-law tails. We show that this results from a stochastic 
cascade of amplification of fluctuations due to transient periods of instability. The exponent of the 
power-law is interpreted as a negative fractal dimension. 
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Intermittency refers to the alternation of long periods 
of rest, followed by bursts of strong activity. This phe¬ 
nomenon plays a crucial role in the dynamics of a wide 
range of physical systems, for example in fluid dynam¬ 


ics 


. dl: 0; astrophysics Q, condensed matter , as well 
as in physiology [3. Theoretical explanations proposed 
in the case of deterministic dynamical systems with a few 
degrees of freedom assume that the system is close to the 
transition to chaotic behaviour . The description of 
intermittency in extended dynamical systems remains a 
challenging issue Sim. 

The importance of noise in the dynamics leading to 
intermittency has been established in pipe flows [l^ and 
spatially extended systems which are close to the insta¬ 
bility onset [HQ, and is likely to play an important 
role in other physical systems as well. This article is 
concerned with the influence of noise on intermittency in 
simple model systems. 

Here, we consider a universal form of intermittency 
arising from the addition of noise to a dynamical system 
which has stable dynamics. It arises for dynamical sys¬ 
tems which are (i) non-autonomous, and (ii) stable, in 
the sense that in the absence of noise, solutions converge 
to a point attractor. It is the interaction between the 
noise and the fluctuating environment which generates 
intermittency, via a mechanism of stochastic amplifica¬ 
tion. In one spatial dimension we consider equations of 
motion in the form: 


X = v{x, t) -\- VWr]{t), (1) 

where v{x,t) is a velocity field fluctuating in space and 
time, ri{t) is a white noise signal with {rjft)) = 0 and 
{il{t)r]{t')) = S{t — t'), and D is the diffusion coefficient of 
the corresponding Brownian motion, which we assume to 
be small {{X) denotes the expectation value of X). In the 
deterministic case {D = 0) the system is characterised 
by its Lyapunov exponent A. If A < 0, the system is 
stable, so all the trajectories of the system converge to 
an attractor [l^ when D = 0. If the system is unstable, 
A > 0, there is typically a strange attractor [^. 


In the presence of noise {D ^0), however, the trajec¬ 
tories do not remain on the attractor of a stable system. 
As a result, the separation of two trajectories, denoted 
as Xx, can have large excursions away from 0, leading to 
tails of the probability distribution function (PDF) Pax 
(T hroughout we denote by Px the PDF of the quantity 
A). 

The PDF Pax may be expected to be well approxi¬ 
mated by a Gaussian distribution. In the case of a stable 
autonomous system the motion in the vicinity of the at¬ 
tractor satisfies x = Xx -\- V2Dr]{t) with A < 0, which is 
an Ornstein-Uhlenbeck (OU) process [l^. In this case 
the deviations from the attractors do have a Gaussian 
distribution, with a variance P/|A|, so the deviations Aa; 
between two trajectories also have a Gaussian PDF, with 
a variance 2 P /1A |. 

We demonstrate here the existence of a generic class of 
models for which the distribution Pax is non-Gaussian 
and has power-law tails over a range of values of Xx: 

Pax - \Xx\-‘^^+°‘\ ( 2 ) 

when Aa: is large compared to yjD/\X\. Note that when 
a > 0, the asymptotic form ([2]) cannot extend all the 
way to Aa; = 0. The large excursions of Xx{t) described 
by the power-law tails are a manifestation of the phe¬ 
nomenon of intermittency, whereby close-by trajectories 
are separated by the combined effect of the noise and 
fluctuating environment. This is illustrated in Fig. [U 
which shows individual trajectories and their differences 
for a simple model discussed in detail below (see Eqs. m 
and fill) !. Figure [2] shows the corresponding power-law 
distribution of Aa;. 

In this paper we explain the origin of this intermittency 
and show how it can be quantified by analysing an equa¬ 
tion which determines the exponent a. We show that 
there are universal aspects to these fluctuations which 
arise because the mechanism for producing the large fluc¬ 
tuations is independent of the mechanism which seeds 
them. The power-law ([2]) is an emergent property, in 
the sense that it arises for generic equations of motion 
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FIG. 1. (Color online), a Set of trajectories of a model 
of colloidal particles in suspension, Eqs. IIOI). Different tra¬ 
jectories separate and recombine, b Intermittent separation 
of pairs of trajectories Ax{t). The parameters are 7 = 0.05, 
4 = 0.08, A = 0.02, D = 10"®. 


which do not themselves contain non-integer exponents. 
We also observe that statistics of the temporal variation 
of the intermittent signal, such as the waiting times dis¬ 
tribution, follows a power-law behaviour. 

Beyond the point at which the underlying system 
becomes unstable, the system has a strange attractor, 
where phase points cluster on a fractal measure As 
A approaches zero from below, the exponent a in Eq. m 
approaches zero from above. When A > 0, the two-point 
correlation function of the strange attractor, g(Ax), has 
a power-law dependence: 


g{Ax) ~ |Ax| 


D2 — 1 


(3) 


where D 2 is the correlation dimension 17|. Equations 
© and ([3]) have the same structure with the exponents 
related by 


a = —D2 


(4) 


so that normalisable distributions of fluctuations corre¬ 
spond to negative values of 02- Equation ([T]) therefore 
gives a physical meaning to a negative fractal dimension, 
but we should emphasise the difference between the inter¬ 
pretation of the two exponents. Equation ([3]) describes 



FIG. 2. (Color online). The probability distribution Pai for 
the model (1101) . for two values of D {D = 3 x 10"^®, full curve, 
and D — 10“®, dashed curves, showing fits by a power law 
IThe parameters are a 7 = 0.3125, 4 = 0.08, 
A = 0.02 b 7 = 0.375, 4 = 0.08, A = 0.02. 


the pair correlation function for a system with no added 
noise, whereas Eq. m describes the tail of the PDF for 
separations when a small noise signal is added. The for¬ 
mal relation between our exponent a and the correlation 
dimension D 2 allows us to use recent advances in com¬ 
puting the correlation dimension 18 - 2^ to quantify our 
exponent a. An alternative definition of negative fractal 
dimension has been offered in 211. 

We now explain the power-law tails by introducing a 
cascade amplification mechanism. Consider the linearisa¬ 
tion of Eq. to give the separation between two nearby 
trajectories: 


Sx = Z{t)5x + 2'/Dr/{t), (5) 

where 

Z{t) = ^{x{t),t). ( 6 ) 

Note that when D = 0, Z{t) is the logarithmic derivative 
of the separation 6x(t), and that its expectation value is 
the Lyapunov exponent: 

m = ^, A = ^hm - dt' Z{t') = {Z{t)) . (7) 

We can think of Z{t) as being an instantaneous Lya¬ 
punov exponent. In the case of autonomous systems with 
an attractor, the attractor must be a fixed point in phase 
space, and Zit) approaches a constant A < 0 as f —>■ 00 . 
In this case the fluctuations are described by an OU pro¬ 
cess and the distribution Pax is Gaussian. In cases where 
the dynamical system is non-autonomous, Z[t) need not 
approach a constant value. If the external driving is a sta¬ 
tionary stochastic process, Z{t) is a fluctuating quantity 
with stationary statistics. The origin of the power-law 
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tails described by ([2|) is that the fluctuations are ampli¬ 
fied during periods when Z(t) > 0. This noise amplihca- 
tion is independent of the initial amplitude, because the 
fluctuating quantity Z{t) acts multiplicatively in Eq. ([5|). 
This leads to a stochastic cascade amplification process, 
whereby large amplitude fluctuations are built up by a 
succession of periods where Z(t) > 0. The power-law 
tail in the fluctuation distribution arises whenever Z{t) 
is positive for some intervals of time, however short. 

In order to quantify this picture, let us consider the 
dynamics of the fluctuations in a logarithmic variable 

Y = In(Ax) . (8) 

Consider the tail of Pax, where the fluctuations are much 
larger than the driving noise, so that the term \/2Drj{t) in 
([5|) can be neglected. In this limit the equation of motion 
for Y{t) is simply Y = Z{t). Because the fluctuations of 
Z{t) are independent of E, the PDF of Y is expected to 
approach a limit which has translational invariance, up 
to a normalisation factor. This implies that: 

Py ~ exp(—oE) (9) 


correlation time, with statistics given by 


{u{x,t)) = 0, 

(u{x, t)u(x', t')) = exp 


(11a) 

5{t — t'). (11b) 


All of the explicit work shown in this paper is based on 
Eqs. and m- This model has two very different 
random elements. The random velocity field u{x,t) is 
the same for all trajectories, whereas the Brownian noise 
r]{t) has a different realisation for each particle trajectory. 
This model has been extensively analysed for D = 0. 
The Lyapunov exponent, obtained in |l5| . was found to 
be negative for sufficiently large values of 7, with the 
attractor not being a fixed point, but a random walk. 
The separation of trajectories was analysed bv 1241 and 
the correlation dimension was investigated in [l9| for the 
case where D 2 > 0. When D ^ 0, we find that the devia¬ 
tions of the trajectories are found to exhibit a power-law 
tail in their PDF, described by Eq. ([2), see Fig. [2 
Consider the linearisation of Eqs. (II2: 


(5a; = -I- 2'/Dr]{t), 

Si) ='y [S{t)Sx — dv] ^ (12) 


where the constant a depends on the statistics of Z{t) 
[23. The corresponding PDF of Ax [Eq. Q] is then ob¬ 
tained by writing the change of variable dP = FydY = 
PAxdAx. Note that when Y —00, or alternatively, 
Ax —>■ 0, the noise term dominates in Eq. ([5]), so Eq. (ED 
does not apply. This prevents any difficulty with the 
divergence of Py- In fact, the role of the noise term re¬ 
duces to providing a natural cutoff at small separations 
(Ax —t 0), thus preventing potential normalisation prob¬ 
lems. It also follows that the exponent a is independent 
of P, provided that D > 0. 

Next we describe a concrete and physically important 
example of a system that produces fluctuations described 
by ©. This is provided by colloidal particles in a turbu¬ 
lent fluid flow, with velocity field u{x(t), t). The motion 
of small particles suspended in the flow is determined by 
viscous drag, which makes their velocity v{t) relax to that 
of the surrounding fluid. When the particles have a den¬ 
sity which is much higher than that of the fluid in which 
they are dispersed, the viscous drag is proportional to 
the difference between the particle velocity and the fluid 
velocity at its current location (limitations of the model 
are discussed in 22, El). The problem studied in this 
paper is based upon a one-dimensional version of this 
model, which includes Brownian motion of the particles: 


X = v V2Dr]{t), (lOa) 

V = ^[u{x,t) — v], (10b) 


where S{t) = |^(x(t),t) is the velocity gradient at the 
position of a particle, which can be modelled by a white 
noise signal with diffusion coefficient 7): 

s{t) = v^at) , v = ^ ( 13 ) 

where ^(t) is independent of r]{t) but has the same sta¬ 
tistical properties. 

When P = 0, from Eqs. (IT^ and (IT^ we obtain the fol¬ 
lowing equation of motion for Z{t) (previously obtained 
in dll): 


Z = --fZ - Z^ + V^Cit). (14) 


We wish to determine the PDF for Y in the form ED- 
From Eq. (1141) we can construct a Fokker-Planck equation 
for the joint probability density p{Y, Z) of Y and Z, and 
seek a solution in the form p{Y, Z) = e'xp{—aY)pz{Z) for 
consistency with (jOD . Following an approach described in 
dl , this leads to a differential equation for a in the form 

Fpz{Z)+aZpz{Z)=Q (15) 


where P is a differential operator defined by writing 


Fp{Z) 


A 

dZ 


{jZ + Z^)+V 


dZ 


p{Z) . 


(16) 


Because d/dZ is a left-factor of F, any normalisable so¬ 
lution of (TTKI) with a 0 must satisfy: 


where 7 is proportional to the viscosity of the fluid. We 
consider a random velocity field with a vanishingly small 


/ OO 

dZ Zpz(Z) = 0. 

-00 


(17) 
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FIG. 3. (Color online) Exponent a defined by ([2]) as a function 
of the dimensionless parameter e defined by Cl. Note that 
a —>■ 1 as e —>■ 0 and that in the vicinity of the critical point 
Ec where A = 0, we have a ~ K{ec — e)- The dashed line uses 
our theoretical value for the coefficient K and the cross shows 
the critical point location. 


Equation (II3 must be imposed on any approximate so¬ 
lutions of (ITSl) constructed by perturbation theory. We 
remark that the integral in (ITT)) is distinct from the Lya¬ 
punov exponent, because pz(,Z) is a distribution of Z 
which is conditional upon the value of Y. 

In order to facilitate the analysis of (USD, we replace 
Z hy a. scaled variable x, and introduce a dimensionless 
parameter e: 


x{t) = 



e = 


With these definitions, Eq. (HH) is replaced by 


A 

dx 


X + ex + 


A 

dx 


p(x) -I- aexp(x) = 0 . 


(18) 


(19) 


Our numerical results indicate that the exponent a, de¬ 
termined by directly computing Pax, is a function of the 
scaling variable e, as illustrated in Fig. [S] We consider 
two different perturbative approaches to determining a 
as a function of e. The first is to make an expansion 
about £ = 0. This can be done by following the method 
discussed in 18| . The series expansion of a has only one 
non-zero term: a = 1, and all of the coefficients of higher 
powers of e are identically zero [^ . The implication is 
that a has a non-analytic dependence upon e. 

An alternative approach is to make a perturbative ex¬ 
pansion about the critical point where the Lyapunov ex¬ 
ponent changes sign. For the model underlying (fTUl) . this 

occurs at £c ~ 1.33_ We follow an approach used in 

[l^ (see also [i^)- To leading order in £ — £c we obtain 


a = K{ec - e) . (20) 

Using the approach described in [l^ and the coef¬ 
ficient K is expressed in terms of finite-dimensional in¬ 
tegrals, whose numerical evaluation leads to K ^ 0.688. 


FIG. 4. (Golor online) PDF of the waiting times T for dif¬ 
ferent values of the noise intensity D. The dashed line cor¬ 
responds to a data fit to the function Pt = NT~'^ exp {—bT) 
with r = 1.12 ± 0.07, and the solid line corresponds to a 
power law with exponent —1.12. The inlet shows results 
obtained by choosing different values of the threshold cth 
(0.002 < Cth < 0.02) for the case of D = 10“®. The parame¬ 
ters of the simulation are 7 = 0.05, A = 0.02 and A = 0.08. 


This value is found to be in good agreement with the 
results shown in Fig. [3] 

We have shown that intermittency in our model leads 
to power-law behaviour of Pax- We now investigate the 
temporal variation of the signal Ax(t). The intermit¬ 
tency of Ax{t) shown in Fig. [1] can be characterised by 
considering the distribution of waiting time intervals T 
over which Ax(t) remains below a defined threshold, say 
Cth. Figure S] shows the PDF Pt of the waiting times 
T for three different cases of noise intensity, namely 
D = 10“^, 10“®, and 10“®, and for the case where 
£ Ri 1.118. We observe that Pt follows a power-law be¬ 
haviour with an exponential decay at long times which 
can be fit to the function Pt = NT~'^ exp {—bT) with an 
exponent r = 1.12 ± 0.07, independent of D. Note that 
varying D only affects the cut-off value of the exponen¬ 
tial tail, i.e. the value of the parameter b: increasing D 
decreases the range of values of T over which Pt has a 
power-law dependence. The inset of Fig. U shows that 
the choice of the threshold Cth does not affect the expo¬ 
nent of the power-law regime, but only modifies b. It is 
important to remark that none of the mechanisms ex¬ 
plaining the power-law distribution of inter-burst times 
observed in other physical systems (see e.g. |13l. 1271 l28l|) 
is applicable in the present problem. 

To conclude, we have explained and characterised a 
class of intermittent fluctuations in dynamical systems. 
They differ from the usual types of intermittency consid¬ 
ered in low-dimensional systems, in that they arise when 
the equations of motion have additive noise, but the un¬ 
derlying dynamical system is stable (i.e. has a negative 
Lyapunov exponent). We have shown using symmetry 
arguments, and the assumption that the instantaneous 
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Lyapunov exponent has positive fluctuations, that the 
intermittency is characterised by a power-law distribu¬ 
tion of the magnitude of the fluctuations, and we have 
discussed perturbative methods for estimating the expo¬ 
nent a. We have also presented evidence that there are 
power law distributions of other quantities, such as the 
inter-burst times, which is a feature common to many dif¬ 
ferent intermittent systems. The stochastic cascade am¬ 
plification mechanism for producing these fluctuations is 
generic, so that they should be observable in a wide class 
of systems. 
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